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In this paper, we extend the first-order post-Newtonian scheme in multiple systems 
presented by Damour-Soffel-Xu to the second-order contribution to light propagation 
without changing the virtue of the scheme on the linear partial differential equations of 
the potential and vector potential. The spatial components of the metric are extended to 
second order level both in a global coordinates {qij/c 4 ) and a local coordinates (Q a b/c 4 ). 
The equations of qij (or Q a b) are obtained from the field equations. The relationship 
between and Q a b are presented in this paper also. In special case of the solar system 
(isotropic condition is applied (gy = Sijq)), we obtain the solution of q. Finally, a 
further extension of the second-order contributions in the parametrized post-Newtonian 
| formalism is discussed. 
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Recently a series of space astrometric missions are proposed, e.g. LISA(Laser Interfer- 
ometer space Antenna),® GAIA (new generation of HIPPARCOS(HIgh Precision PARal- 
lax Collecting Satellite) for attaining angular accuracy of about 10 -6 arcsec in few hours of 



integration),® SIM(Space Interferometry Mission),® DIVA(Double Interferometer for Visual 
^ ■ Astrometry) , ® FAME (Full-sky Astrometric Mapping Explorer),® ASTROD(Astrodynamical 

Space Test of Relativity using Optical Devices)®, Mini-ASTROD® and so on. LISA and 
GAIA are taken as cornerstones of ESA(European Space Agency). SIM and FAME have 
been planned by NASA(National Aeronautica and Space Administration) and DIVA is in 
the progress in Germany. Especially Mini-ASTROD is proposed by our Chinese colleagues 
and is under Phase A Study.® In all of above missions the light propagation has to be calcu- 
lated to a more precise level, e.g. in Mini-ASTROD the precision is about 10 -13 , which is the 

second order post Newtonian(2PN) level (2PN precision for the sun is 10~ 12 f-^r-J , where 
Rq and r are the radius of the sun and the closest distance from light trajectory to the center 
of the sun). Therefore the second order contribution to light propagation is required to be 
considered. On other hand, when the light pass nearby the sun or any planet, the potential 
produced by 2^ multipole moment is® V mono (l/r) N , where V mono is the potential created 
by the mass monopole, I the length deviated to monopole and r the distance of the action. 
For the quadrupole of the sun (N = 2) (l/r) 2 < e 2 , where e 2 is the small parameter of post- 
Newtonian expansion (nearby the sun e 2 ~ 10~ 6 ), / can be estimated from the oblateness of 
the sun®. Therefore if we consider 2PN problem, the 1PN quadrupole has to be taken into 
account. When we calculate relativistic multipole moments and astrodynamics (equations of 
motion for plannets), we have to use the scheme presented by Damour-Soffel-Xu^ 10-13 ] , in 
which a complete 1PN general relativistic celestial mechanics for N arbitrarily composed and 
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shaped, rotating deformable bodies is described. Their scheme is widely abbreviated as DSX 
schnW 14 ] . In their first paper! 10 ! ; they give out the clear relation between potential w & vector 
potential wi in a global coordinates (they are directly related with the global metric tensor) 
and W & W a in a local coordinates. It is so called "Theory of Reference Systems", which 
can solve the difficulty caused by what the multipole moments is calculated in every local 
coordinates, but the influence of multipole moments on the light-ray trajectory is considered 
in the global coordinates. Their paper II^ 11 ! first obtained the complete and explicit laws of 
the motion for N arbitrarily composed and shaped bodies in the form of an infinite series 
which make the relativistic correction of multipole moments on astrodynamics (Ephemeris) 
to be possible. Therefore DSX scheme is the best framework for taking into account the 1PN 
multipole moments. However DSX scheme is the first post-Newtonian (1PN) approximation 
for particle motion. The metric tensor in a global coordinates is written in the form 



<m - -cx P (-^)+o(6), (l) 



c 



Aid- 

90i = -^ + 0(5), (2) 
9ij = <%exp(|^) +0(4), (3) 

and they satisfy 

9oo9ij = Sij + 0(4) , (4) 

where 0(n) means 0(c~ n ). Eq.(4) is the conformal isotropic condition (see Eq.(4.5) of 
Ref.flO]). Potential w and vector potential Wi satisfy linear Partial Differential Equations(PDE). 
The metric in a local coordinates has a similar form as in the global coordinates, one only 
needs to change the small letters into Capital letters. It has been pointed out already, that 
DSX scheme has to be extended to 2PN, then it is possible to discuss 2PN contribution' 15 ] 
to light propagation, which means the Eq.(3) to be extended to 0(6). 2PN contribution for 
lightray has been discussed for a long time (early 80's) by Epstein and Shapiro (1980)' 16 ], 
Richter and Matzner(1982)[ 17 ] , and by others (later). But all of them consider only in one 
global coordinate system, therefore they can not calculate the relativistic contribution from 
multipole moments which should be calculated in the local coordinates. As we know the 
relativistic theory of reference system established only after 199l[ 10 ], although it is only 1PN 
approximation. 

In this paper we will first extend the metric g^ (and G a b) in DSX scheme to 0(6), and 
deduce the equations satisfied by w, Wi and q^ (also for W, W a and Q a t, in every local 
coordinates) from Einstein field equations, where we introduce new quantities (and Q a b in 
the local coordinates) for metric gij. Then the relationship between potentials in the global 
and the local coordinates is deduced also. After that, the special case of the solar system 
(Qij = <%<?)' 17 ' is discussed, and an integrable solution of q is presented. Finally we give a 
brief discussion. 

All of symbols in this paper and conventions are the same as in DSX schemel 10 " 13 ]. To 
apply DSX scheme to discuss the second order contribution to light propagation, g^ (or G a b) 
should be extended to 0(6) as we mentioned before. 

The metric tensor in a local coordinate system now take the form 

Goo = -exp(-— ) +0(6), (5) 



c- 



G 0a = -^ + 0(5), (6) 
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G ab = 6 ab exp(™) +% + 0( 6 ). (7) 



The contravariant metric reads 



\ c 2 J c 4 



(2W\ 



G 00 = -exp^j+0(6), (8) 

G 0a = _l^ + 0(5)) (9) 
c? 

G ab = 5 a6 exp(-^)-^ + 0(6). (10) 

The metric tensor in the global coordinate system has a similar form, just change G a p, W, 
W a and Q ab by g^, w, Wi and 

The relations between w , W{ , ffoj and , W a , Qab will be studied in the following. 
In DSX scheme, since the metric in the global coordinate system as well as in every local 
coordinate systems has the similar form, so that, if we deduce any equation in one local coor- 
dinate system, then we have it in every local coordinates as well as in the global coordinates. 
Therefore we need only deduce the field equation in a local coordinate system. The spatial 
conformal isotropic condition Eq.@ is revised as 

G ab G 00 = -5ab - % + 0(6) . (11) 

If we attribute Qab/c 4 to 0(4), it is just the form in DSX scheme. Therefore Ea.pi|) is an 
extension of Eq.ljlj). Q ab also can be taken as a spatial anisotropic contribution in second 
order. 

The metric tensor goi is taken as e 3 here, in fact it is only as small as e 4 because of 
the relatively slow rotation rate in the sun I 17 ' 18 !. To keep the uniformity with DSX scheme 
therefore we take such formula of goi as well as for Qo a - 

The Christoffel symbols can be calculated from the metric tensor 

r oo = -^ + 0(5), 

rl = -^ + o(6), 

r 00 = -— + ^4 ^ + °( 6 )> 

T° ab = 5 ab ^ + ^W {a , b) + 0(5), 
r o& = -^M + ^ afe + 0(5), 

n c = ^abW, c + 5 ac W^-5 bc W, a ) + ^(Q abtC + Q ac ^-Q bCta ) + 0(6), 

where two indices enclosed in a parentheses is denoted as symmetrization, and in a square 
bracket means antisymmetrization. 

Then we can deduce the Ricci tensor 

^oo = _Y^K_}_( 3dtdtW + 4dtddWd) + o(6), (12) 



R 0a = -^(v 2 W a -d a d d Wd-d t d a w)+0(5), (13) 

45 ab W V 2 W + S ab W >tt + 4W (aib)>t - 2W, a W, b 

+ 0(6), (14) 



2, 



R ab = -^5 ab V 2 W+ 1 



c 2 a ° c 4 



+ -^{Qad,bd + Qbd,ad ~ Qab,dd ~ Qdd,ab) 
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and the scalar curvature 
R = -\v 2 W + ^-(4WV 2 W + 6W, tt + 8W a>at - 2W, a W, a + Q ab , ab - Q aa , bb ) + 0(6) . (15) 
Prom Einstein field equation we can derive equations of W, W a and Q ab 

V 2 W + \{3W,u + Ad t d a W a ) = -4vrGS + 0(4) , (16) 

V 2 W a - d a 8 b W b - 8 t d a W = -47rG£ a + 0(2) , (17) 
-2*06 + 4(W (aib)>t - <U^ (cM)jt ) - 2W a W )6 + foW^d 

+ ^(Qad,bd + Qfed,ad ~ Qa6,dd _ Qdd,ab ~ KbQdc,dc + 5abQdd,cc) = 8irGT ab , (18) 

where £ = (T 00 + T aa )/c 2 and S a = T 0a /c (the same definition in DSX scheme). Eq.^ and 
(|17|) are linear PDE, which can be solved in certain suitable gauge conditions, then Eq. (|18[) 
will be solved also. We shall point out that, T ab can be expressed by E and S a (to see 
Eq.(A13) of Ref.[19]). Therefore it is self-consistent within the framework of DSX scheme. 
Similar field equations in the global coordinate system can be obtained easily by substituting 
capital letters by small letters. 

From Eqs.(16)-(17) we see that Q ab does not appear in the field equations of W and W a , 
and they are the same as ones in DSX scheme (keeping linear equations). The solutions of 
W and W a related to relativistic multipole moments are still valid as before in DSX scheme. 

Now, we shall discuss the transformation relations of the potential and the vector potential 
from the global coordinate system to a local coordinate system and vice versa. The coordinate 
transformation law reads 

n V dx^ dx u p 
9 dX<*dX? ' 1 > 

Considering g 00 , we obtain 

/ 1V 2 \ A r 2 

w = [ 1 + J W + ^V a W a + - ln(4)4) - A °a A a) + 0(4) , (20) 

where A£ = dx^/dX a , V a = Rfv 1 (Rf is a slowly changing rotation matrix and v l is the 
coordinate three- velocity of the central world line measured in the global coordinate system) . 
For g° l , it turns out to be 



c 
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Wl = v *w + KW a + -(A» Al - AIA\) . (21) 

In Eqs.(|2()jl and (|2*T|l . Qij does not enter these equations, and the relations between w, Wi 
and W, W a are the same as in DSX scheme (to keep linear relation). As we know, all of the 
discussion about the theory of reference systems are based on the linear PDE of potentials 
(Ea. ()16j) . 1)17(1 ) and linear relation between potentials in the global coordinates and in the 
local coordinates ((|2()jl and (|21j) ). Since these equations are the same as before (in DSX 
scheme), the theory of reference system in DSX scheme (for potential and vector potential) 
is valid in this paper also. 
Finally, for g lJ , we get 

qij = -2W(2V 2 5ij - RiRlv a V b ) - 8^W a <% + &R%B$W a V h + R a RlQ ab 

+2Wc 2 (AiAi - 5 t3 ) + c 4 {44 - AlAi + <%[1 - ln(A A° - A° a A° a ))} . (22) 

Here we should emphasize that Eq.(|22jl has only formally given for completeness. When we 
consider the problem of the light propagation in the solar system Ea. (|22() is unnecessary to 
be considered which we will explain later. 
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To discuss the gravitational influence on light propag ation for Mini-ASTROD^ (1(T 13 ) 
we consider the solution in the solar system now. Potential W and vector potential W a can 
be obtained in the same as in Ref.[10] by means of B-D multipole moments' 20 ' in every local 
coordinates. Therefore we shall not repeat it again. Q a b is rather complex. Fortunately for 
the problem of the light propagation in the solar system it becomes much simple. Since the 
velocity of the relative motion is small inside the sun (f 2 /c 2 < e 2 ), especially the shape of 
the sun is close to a monopole, therefore if the origin of the coordinate system is taken as 
the center of the solar mass (dipole always equal to zero), the quadrupole terms (and higher 
multipole moments) in 2PN contribution on Q a b can be ignored, i.e. Q a b = S a bQ, as well as 
in the global coordinates of the sun 

qij = Sijq. (23) 

As for other local coordinates of every planets 2PN contribution is too small (10~ 19 — 1CP 18 ) 
to be calculated. 

Substituting Eq.JH into En. fTCty and taking PN gauge (3d t W + 4d a W a = 0(2)), Q 
satisfies the equation as following: 

Q + \w 2 \ = 8irG (T aa - ^) . (24) 
1 J ,dd V 2 / 

The right hand side of above equation has a compact support and the solution of Eq. (|24[) is 



Q = -/ 2G(T ° J -^ /2 W-V. (25) 



From Ea .1)25 [) it is clear that Q itself does not have compact support neither is linear in W 
(because of the last term). 

In the global coordinates of the solar system, we have a similar equation and solution. 
Since there is very small difference between the global coordinates of the solar system and the 
local coordinates of the sun (additional 10 -3 in c -4 level) we can always ignore the distinction 
between q and Q of the sun. This is the reason why Eq.(|22|) is unnecessary to be considered. 

We have extended DSX scheme to 2PN level approximation. In the solar system we obtain 
the solution of q and Q. Our extension keep the main virtue of DSX scheme (linear equations 
of the potential and the vector potential and linear relationship between w, w l and W, W a ). 
Here we discuss further on our results. 

This extension may also be used in parametrized post-Newtonian PPN formalism, which 
would be found a widely use for the relativistic experiments and observation! 21 ]. For DSX 
scheme (1PN) the PPN extension is written as' 22 ' 

5 oo = -1 + — - ~Pw 2 + 0(6) , (26) 
c z c 4 

2(1+7) • 

9ij = (% (l + + 0(4) . (28) 

Our extension for parametrized second post Newtonian (PP 2 N) formalism for the solar system 
is taking the form 

On, 9 

(29) 
(30) 

2e f GT SS , 3 f 




c 4 



/^-^x + l/^dx 3 +0(6), (31) 
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where T is the energy-momentum tensor in the global coordinates. When j = f3 = 5 = e = 
£ = 1, then the metric tensor return to the case in general relativity. If we ignore all of 1/c 4 
terms in g^j, our results agree with the one of Klioner and Soffel J 22 l 

We would like to thank Prof. Ni Wei-tou and Prof. Tao Jinhe for the useful discussion. 
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